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Abstract. Electron tunneling through small metallic islands with low ca- 
pacitance is studied. The large charging energy in these systems is respon- 
sible for nonperturbative Coulomb blockade effects. We further consider 
the effect of electron interactions in the electrodes. In junctions with high 
resistance compared to the quantum resistance transport can be described 
by sequential tunneling. If the resistance is lower, quantum fluctuations, 
higher order coherent processes, and eventually resonant tunneling become 
important. We present a path integral real-time approach, which allows 
a systematic diagrammatic classification of these processes. An important 
process is "inelastic resonant tunneling", where different electrons tunnel 
coherently between the electrodes and the island. Physical quantities like 
the current and the average charge on the island can be deduced. We find a 
strong renormalization of the system parameters and, in addition, a finite 
lifetime broadening. It results in a pronounced broadening and smearing 
of the Coulomb oscillations of the conductance. These effects are impor- 
tant in an experimentally accessible range of temperatures. The electron 
interaction in the electrodes is modeled by a Luttinger liquid. It leads to 
non-analytic kernels in the effective action. The diagrammatic expansions 
can be performed also in this case, resulting in power-law current-voltage 
characteristics. 
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1. Introduction 

Electron transport through mesoscopic metahic islands coupled to elec- 
trodes has been the subject of extensive research [1-3]. The small size and 
capacitance of these systems imply a strong Coulomb interaction energy 
Ech, which gives rise to a variety of single-electron phenomena. At low tem- 
peratures tunneling can be suppressed (Coulomb blockade). The charging 
energy can be tuned by gate voltages. When the energy difference between 
adjacent charge states Ech{n ± 1) — Ech{n) is lower than the temperature 
kT or the bias voltage eV, a current flows through the system. As a conse- 
quence, the conductance shows a series of peaks as a function of the gate 
voltage Vg (linear response) and further structure at larger bias voltage V 
(nonlinear response). The classical description is sufficient as long as the 
resistance Rt of a single barrier is much higher than the quantum resis- 
tance Rk = h/e^, i.e. for ao <C 1, where ao = Rk / {^t^^ Rt) ■ In this regime, 
transport occurs in sequences of uncorrelated tunneling processes. The rates 
can be obtained in lowest order perturbation theory in the tunneling am- 
plitudes. They enter a master equation [4-8], from which the probabilities 
for different charge states n and the currents can be calculated. 

For very low temperatures or when the dimensionless conductance ao 
is not small the classical description breaks down. Quantum fluctuations 
and higher order coherent tunneling processes become important [9-15]. 
This includes cotunneling, where two electrons tunnel coherently in differ- 
ent junctions, thus avoiding the Coulomb blockade. Furthermore, resonant 
tunneling, where electrons tunnel coherently back and forth between the 
island and the electrodes, play a role. A description is called for, which 
allows a systematic classification of all these processes. In comparison to 
the well known phenomenon of resonant tunneling of single electrons one 
encounters here two complications. One lies in the fact that the metallic 
system contains many electrons. With overwhelming probability different 
electron states are involved in the different transitions of the coherent pro- 
cess. The second arises since the Coulomb interaction is strong and, hence, 
cannot be accounted for in perturbation theory. 

In the present article we develop a systematic diagrammatic technique 
to identify the processes of sequential tunneling, inelastic cotunneling and 
resonant tunneling. We study the time evolution of the density matrix. In 
an earlier paper [16] we have formulated the problem, after a separation of 
charge and fermionic degrees of freedom, in a many-body expansion tech- 
nique. Here we reformulate it in a real-time path-integral representation. 
The latter is well known from the studies of dissipation in quantum mechan- 
ics. Caldeira and Leggett [18], following methods pioneered by Feynman and 
Vernon [17], studied the problem of Ohmic dissipation. Dissipation associ- 
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ated with tunneling of electrons has been investigated in Refs. [19, 5]. An es- 
sential step in the present work is a transformation of the Feynman- Vernon 
functional for electron tunneling from a phase (coordinate) representation 
to a charge (momentum) representation [5]. The path- integral method is 
well-suited to account for the strong correlations due to the Coulomb in- 
teraction in a perturbative or nonperturbative analysis in the tTinncling. 
In contrast, usual Green's function techniques cannot be used since the 
Coulomb interaction has to be retained in the unperturbed Hamiltonian. 
The same problem arises in the context of local, strongly correlated Fermi 
systems like the Kondo and Anderson model [20, 21]. For these systems dia- 
grammatic techniques have been derived by Barnes [22] from a slave-boson 
description. Rammer [23] developed a graphical density-matrix description 
for the dynamics of a particle coupled to a heat bath. The problem can also 
be formulated in terms of Liouville operators [24] . 

As examples we will study the electron box and single electron tran- 
sistors. The box consists of a metallic island coupled via a tunnel junction 
with capacitance Cj to an electrode. It is further coupled capacitively (Cg) 
to a voltage source Vq- The charging energy is Ech{n) = {ne — Qg)'^/ (2C), 
where C = Cj + Cq denotes the total capacitance of the island. It depends 
on the number of excess electrons n on the island and on the continuously 
varying external charge Qg = CgVg- The single electron transistor consists 
of a metallic island which is coupled by two tunnel junctions to two elec- 
trodes (see Fig. 1). Here a transport voltage V = Vl — Vr drives a current. 
The island is further coupled capacitively to a gate voltage Vq- The charg- 
ing energy of this system depends again on the number of electrons n on the 
island, Ech{n) = {ne-QG f /{2C). Here C = Cl + Cr+Cg is the sum of the 
two junctions and the gate capacitances and Qg = CgVg + ClVl + CrVr. 
In tunneling processes we further have to account for the work done by the 
voltage source (see below). In both examples the total island capacitance 
defines the scale for the charging energy Ec = e^/2C. 

The tunneling is described by tunneling Hamiltonians. We consider 
"wide" metallic junctions, which implies that there are many transverse 
channels. As a result "inelastic" higher order tunneling processes, involving 
different electron states for each step, dominate over those higher order pro- 
cesses which involve the same state repeatedly. Accordingly, in the effective 
action description presented below only simple loop diagrams arc retained. 

At low temperature the electron number uIQg) in the electron box 
takes the value which minimizes the charging energy. It increases in steps 
of unity when Qg is increased. Tunneling processes in the transistor are 
only possible in lowest order perturbation theory if the electrochemical 
potential of one electrode is high enough to allow one electron to enter the 
island, say eVi > Ech{n+l) — Ech{n) , while the electrochemical potential of 
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Figure 1. Equivalent circuit for the SET-transistor. 

the other electrode allows the next tunneling process to that electrode, i.e. 
Echin + V) — Ech{n) > cVr. Within the window set by these two conditions 
the current is ARtI = V - A{Qg - (n + l/2)e]V(CV). This implies that 
the conductance G{Qg) = dl/dV shows, as a function of Qg, an e-periodic 
series of structures of width CV with, at T = 0, vertical steps at its limits. 
At finite temperature the average value {h^Qg)) and the steps in the linear 
and nonlinear conductance are washed out. 

Quantum fluctuations further wash out these steps, even at zero tem- 
perature. We will describe these processes diagrammatically. In the most 
interesting case we can re-sum the diagrams and obtain closed expressions 
for the stationary density matrix and the spectral density describing the 
charge excitations of the system. Our main findings are: (i) After a resum- 
mation of the leading logarithmic terms in In (E^c/ja;!) we find a renor- 
malization of the energy and dimensionless conductance, (ii) The coherent 
processes where electrons tunnel an arbitrary number of times between the 
leads and the island (resonant tunneling) give rise to a broadening of the 
charge state levels. Both effects need to be retained in a conserving theory 
which obeys sum rules and current conservation. Both effects are observable 
in a real experiment at accessible temperatures. 

In order to study the effect of electron interactions we consider in the 
last section a model where the electrodes are assumed to be Luttinger 
liquids [25] . Due to the interactions the single particle density of states has 
a power-law asymptotic at low energies, and Fermi liquid theory does not 
apply anymore. For a quantum wire with an arbitrarily small barrier this 
leads to a suppression of transport at low energies [26, 27]. We find that 
the effective action in the Luttinger liquid - normal metal junction has the 
same form as in the metallic case. However, the kernels are modified, taking 
a power law form with exponents depending on the interaction strength 
in the Luttinger liquid. The diagrammatic expansions can be performed 
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and re-summed also in this case. The resulting /-^-characteristics show a 
power-law dependence. 

2. Real time evolution of the density matrix 

The description of the single electron transistor is based on the Hamiltonian 
H = HL + HR + Hi + H,h + Ht. (1) 

Here Hl = J2k ai^k + ^VLjf^ka^ka describes noninteracting electrons in the 
left lead, with similar expressions for the island (with states denoted by 
q) and the right lead. The electrodes are treated as reservoirs, i.e. they re- 
main in thermal equilibrium. The Coulomb interaction Hch is assumed to 
depend only on the total charge on the island, as expressed by the charg- 
ing energy Ech{n) introduced above. Charge transfer is described by the 
standard tunneling Hamiltonian, e.g. the tunneling in the left junction by 

Ht,L = '^kqcl^Cqa + h.C. . (2) 

k,q,cr 

The tunnel matrix elements T^g are considered independent of the states 
k and q. They can be related to the tunnel restistances Rl/r of the junc- 
tions, 1/Rl/r = (47re2/n)iVi/^(0)iV/(0)|Ti/^P, where Nr{0) arc the den- 
sities of states of the the electrodes and the island, r = L. R, I. We further 
assume that the junctions, although small, still accommodate many trans- 
verse channels. From a comparison of Andreev reflection and single electron 
tunneling in small normal-superconducting junctions we have concluded 
that the number of channels is N^h > 10^. 

The quantum mechanical many-body problem of electrons coupled by 
Coulomb interactions can be reformulated in a path integral representation. 
In order to handle the interaction one performs a Hubbard-Stratonovich 
transformation. It introduces collective variables ipr, which are the quan- 
tum mechanical conjugates of the charges in the electrodes or island. The 
capacitive interaction between electrons is replaced in this way by an in- 
teraction of the electrons with the fields Vr = htpr/e. After this stage the 
electronic degrees of freedom can be traced out. The next step involves 
an expansion in the electron propagators. Since we consider wide jimctions 
only simple loops need to be retained [28] . Their iteration introduces in each 
order a factor N^h, hence they dominate over more complicated higher or- 
der loops. After the trace has been performed the system can no longer be 
described by a Hamiltonian. Instead we deal with a reduced density ma- 
trix /o({<^r,i}) Wr,2})- It depends on an effective action, which is expressed 
in the phases ipr,cr corresponding to the forward and backward propagator 
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£7 = 1, 2. The structure of the theory is famihar from Refs. [17, 18], where a 
quantum system coupled to a harmonic oscillator bath has been considered. 
The analogous model describing the effect of electron tunneling has been 
presented in Refs. [19, 5]. We, therefore, do not present the derivation of 
the steps here, rather we quote only the result. 

For transparency we describe the formalism for a single junction. The 
time evolution of its reduced density matrix is given hy {ti = k = 1) 

r ffif rve{ 

p{tf;'Pif,^2f)= dipudip2i P(/9oo(U) / 'Dip^{U) 
J J(pii Jfe) 

(3) 

Here ipi and (p2 refer to the forward and backward time evolution. The 
effective action is given by 

S[^Pl, ^2] = Sch[<Pl] - Sch[^2] + St[ipi, ip2] , (4) 

The first term represents the charging energy Schif] = /// dt^ ((^/e)^. Elec- 
tron tunneling is described by [19, 5] 

St[ipi, ip2] = 47rz Yl r dt f - t') cos[^a{t) - ^^t')] , (5) 

ct,ct'=1,2 

where a'^'^ are given in Fourier space by 

a'''\u) = {-iy+^a-{u) , a^'^iu) = {-lfa+{u) (6) 
and a^(uj) = iaocv ; — — . (7) 

The tunneling term couples the forward and backward propagators. This 
arises in the step where the microscopic degrees of freedom are eliminated. 

An important step for a systematic description of tunneling processes 
is the change from the phase to a charge representation, accomplished by 

/9(t/;ni/,n2/) = ^/o(ti;ni,ni) / d(pifd(p2fd(piid(p2i 

rfif rfe{ r r 

/ Vip^(U) / P^e(U) / P\oo(U) / P\g(U) 

X exp (^-iriiipii + inif(pif - iSch[ni] + i J dtn\ipi 
X exp (^iniip2i — in2fip2f + iSch[n2] — i J dtn2(p2 
X exp{47ri ^ dt t dt' a"^"' {t - t') cos[<^<,(i) - ^Pa'{t')]} . (8) 
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Here Schin] = J^^ dt-^{ne)^. In systems with discrete charges the integra- 
tions over If include a summation over winding numbers [5]. We further 
assume that the initial density matrix is diagonal in the charges. In the ab- 
sence of tunneling the bare forward and backward propagators involve only 
the charging energy exp(zbi5c/i[n]). We expand the tunneling part of the 
action exp{iSt) and integrate over ip. Each of the exponentials exp[±iipfj{t)] 
causes a change of the number of electrons on the island by ±1 at time t on 
the forward or backward branch, a = 1 or 2, respectively. These changes 
occur in pairs and are connected by a line representing a'^'"' {t — t') . The two 
correlated transitions can occur on the same branch cr = cr', or on different 
branches. 

The latter are of particular interest. Imagine we started in a state with 
n charges p{ti) = \n){n\. Then each transition, described by exp{i[ipi{t) — 
f2it')]) changes the charge on both branches by +e and the density matrix 
acquires a finite value also for states |(n -|- l)e). After integrating over the 
two times t and t', limited hy ti < t' < t < tf, we find 

(n + l\p{tf)\n + 1) = {tf - U)27ra+{AE^h{n)) , (9) 

where AE^hin) = E(.fi{n + 1) — Ef.fi{n). Obviously we can interpret the 
coefficient of the time difference as transition rate, and indeed we reproduce 
the well-known single electron tunneling rate. 

In the following we consider the single electron transistor. Each capacity 
Cr introduces a charging term Sch, where r = L, R,G denotes the reservoirs 
on the left, right and the gate, and each tunnel junction introduces a term 
St, all depending on the appropriate phase difference ipi — ipr- The phases in 
the reservoirs r are assumed to be well-defined quantities cpr = eVrt without 
fluctuations. The contribution of the charging energy which is related to the 
work done by the transport voltage can be accounted for in the tunneling 
lines by af{ijO — eV^), and the charging energy becomes Ech{n) = (ne — 
Qg)V{2C). 

Each term of the expansion can be visualized by a diagram. In Fig. (2) 
we show several important processes. There is a closed time-path consist- 
ing of two horizontal lines joined at tf. They correspond to the forward 
propagator from ti to tf (upper line) and the backward propagator from t / 
to ti (lower line). Along the time-path we arrange vertices. They are con- 
nected in pairs by (dashed) tunneling lines, cither within one propagator or 
between the two propagators. We draw all topological different diagrams 
with directed tunneling line and evaluate the diagrams according to the 
following rules: 

1. Assign charge states n and the corresponding charging energy to each 
element of the propagators. Elements of the forward (backward) prop- 
agator between t and t' < t carry factors e-xp[^iEch{n){t — t')]. 
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Figure 2. Example of a diagram showing various tunneling processes: on the left se- 
quential tunneling in the left and right junctions, then a term which preserves the norm, 
next a cotunneling process, and on the right resonant tunneling in the left junction. 

2. Each vertex represents an exponential exp[ibz(^o-(0] tunneling 
contribution to the action. It changes the charge from n to n it 1. 

3. Pairs of vertices are connected by a directed tunneling line a^{t — 
t') [a~ {t — t')] for the electrodes r = L, i?, if the line of is running 
backward [forward] with respect to the closed time-path. The charge 
increases (decreases) along the time-path by 1 if a tunneling line comes 
in (goes out). 

4. Each diagram carries an prefactor (— i)^( — I)'", where M is the total 
number of vertices and m their number on the backward propagator. 

5. Integrate over the internal times and sum over the reservoirs. 

In order to calculate stationary transport properties it is convenient 
to change to an energy representation. This is achieved by the following 
transformation. In each diagram we order the times from left to right and 
label them by tj, irrespective on which branch they are. We further set 
ti = —oo and = 0. We then encounter integrals of the type 

f dh f dt2... f dtM e^*ie-^^-^i(*2-*i)e-^^-^2(*3-t2) . . . ^-iAEM{-tM) 

J —OO Jtl JtM-1 

= 1 1 1 

AEi — if] AE2 — irj AEm — iv 

Here AEj is the difference of the energies of the upper and lower propagator 
and - if present - the frequency of the tunneling line within the segment 
limited by tj and tj+i- The convergence factor e^^^ {rj 0"*") is related 
to an adiabatic switching of the tunneling term Hf. The rules in energy 
representation read: 

1. Draw all topological different diagrams. These are the same as in time 
space. In addition to the charging energy assigned to the propagators 

we assign a frequency lo to each tunneling line. 

2. For each segment derived from tj < t < tj+i with j > 1 we assign 
a resolvent ^ where AEj is the difference of the energies of the 
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forward and backward propagator, plus the sum of the frequencies of 
the tunneling lines in the given segment. The latter have to be taken 
positive for lines from the left to the right and negative for lines from 
the right to the left. 

3. The pref actor is given by (—1)"*"'"', where m is the total number of ver- 
tices on the backward propagator and I the total number of resolvents. 

4. For each coupling of vertices we write a^{uj) [a~ (uj)], if the tunneling 
line of reservoir r is going backward (forward) with respect to the 
closed time-path. 

5. Integrate over the frequencies of tunneling lines and sum over the reser- 
voirs. 

3. Master equation and stationary probabilities 

We denote the sum of all diagrams by 11^^'",^, where the indices indicate 
the left and right charge states on the two branches, a = 1,2. They can be 
expressed by an irreducible self-energy part S"^'^,^, defined as the sum of all 

diagrams where any vertical line cutting through them crosses at least one 
tunneling line. The time evolution of the density matrix can be expressed 
as an iteration in the style of a Dyson equation, illustrated in Fig. (3) 

^nun[ ^ Yi(^)2s^ A + y n"^'";^ eI/I nW"? . (10) 

Here H^^"^ is the free propagator. 

We start from an density matrix which is diagonal, p{ti,ni,n2) = 
P^^\ni)Sni,n2- This means that the density matrix remains diagonal except 
during transitions described by S. Hence, we consider S„ „' = E^'", . We 
identify the solution of Eq. (10) - after multiplication with P^^\n) from 
the left - as the stationary distribution function J2n P^^'' i''^)^n'n' ~ P^^ij^')- 



Figure 3. The iteration of processes for II, describing the time evolution of the density 
matrix. 
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Hence Eq. (10) reduces to 

P^\n') = P(0) (nO + Y: P^*(n")Sn",n'ni?„, . (11) 

n" 

The last term in Eq. (11), n^^^, describes a propagation in a diagonal state 
(i.e. equal energies on the forward and backward propagator) which does 
not contain a tunneling line. According to the rules in the energy repre- 
sentation this introduces a factor 1/ir], and in the limit r] — > 0+ we find 
X]n' P^^i^')^n':n = . By attaching the rightmost vertex of each diagram 
S to the upper and lower propagator we can show that J^n' ^n,n' = , 
which allows us to rewrite Eq.(3) also in the form 

= -P'\n) T.n,n' + E ^'*(«')Sn',n • (12) 

We recover the structure of a stationary master equation with transition 
rates given by Sn'.n- In general the irreducible self-energy E yields the rate 
of all possible correlated tunneling processes. We, furthermore, see that the 
stationary solution P^^{n) is independent of the initial distribution P^^\n). 

For illustration we evaluate now all diagrams which contain no over- 
lapping tunneling lines. After each transition the system is in a diagonal 
state of the density matrix and a probability can be assigned. Examples 
are visualized on the left hand side of Fig. (2). This means we include 
those processes which are also described by the master equation with rates 
obtained in lowest order perturbation theory. In the limit considered the 
irreducible part S*^^) consists of two classes of diagrams, those where a tun- 
neling line connects two vertices on one propagator (which implies that 
n[ = rii) and those where it connects the two propagators (n^ = ± 1, 
with ni — 712 = n'l — 713) (see Fig. (4)). According to the rules the rates are 

sS±i = 2vria±(±A£;±) , E« = -27ri ^ "^(±A£;± ) (13) 

± 

where AE^ = E^h{n ± 1) - E^h{n). 

4. Higher order tunneling processes 
4.1. COTUNNELING 

In situations where single electron tunneling is suppressed by Coulomb 
blockade the lowest order contribution to the current arises due to cotun- 
neling. It is described by a diagram with tunneling processes in the left 
and in the right junction, where the corresponding lines aiiti — t'^) and 
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Figure 4- Representation of the self-energy E'^', defined to contain no overlapping 
tunneling lines. Only one representative of each class is shown, the remaining ones are 
obtained by changing the direction of the arrows and exchanging the position on the 
forward and backward propagator. 

0(R{tR~t'fl) overlap in time. This means there is no intermediate time where 
the density matrix is diagonal, which would describe sequential tunneling. 
An example is shown in Fig. (2). Performing the integrations we find for 
the rate 



plus a divergent term. This divergence is canceled if we take into account 
consistently all processes of the same order. This includes a process where 
the island has already the charge (n zt l)e and one tunneling event brings 
it back to ne. The corresponding terms are products of a single electron 
tunneling rate of one junction and the stationary probability Pn±i in hrst 
order in ao- Adding all these terms we obtain the result given above. 

4.2. RESONANT TUNNELING 

The perturbative approach breaks down at low temperatures or for large 
values for the dimensionless conductance ckq. Specifically we will show that 
the classical master equation is valid only for In (^^) ^ 1, whereas for 
larger values resonant tunneling processes become important. 

To proceed we have to find a systematic criterion which diagrams should 
be retained and summed. For this wc note that during a tunneling process 
the reservoirs contain an electron excitation. Two parallel tunneling lines 
imply two such excitations. Our criterion is that we take into account only 
those matrix elements of the total density matrix, i.e. reservoirs plus charge 
states, which differ at most by two excitations in the leads or (equivalently) 
in the island. Graphically this means that any vertical line will cut at most 
two tunneling lines. 




(14) 
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We will concentrate here on situations where only two charge states 
with n = 0, 1 need to be considered. This is sufficient when the energy 
difference of the two charge states Ag = E^hi'^^) — EchiO) and the bias 
voltage eV = cVl — cVr are small compared to Ec, which is the energy 
associated with transitions to higher states. The combination of the two 
restrictions, diagrams with at most two parallel tunneling lines in the two- 
state problem, implies that no diagram contains crossing tunneling lines. 
In this case we can evaluate the irreducible self-energy analytically. 

In order to sum all diagrams which may contain up to two parallel 

tunneling lines we introduce a diagram labeled by ^^j'"' (^j^^) (see Fig.(5)). 
It has an open tunneling line associated with tunneling in the junction r 



Figure 5. a) Self-consistent equation for 0^tj). A summation over the electrodes r and 
the direction of the tunneling lines is implied, b) Representation of the self-energy E„,i 
within our approximation. 

carrying the frequency uj. Consequently it remains in an off-diagonal state 
at one side. For the two state problem we need only 

<Pni^) = Kiir,uj) (15) 

with n = 0, 1, for which we can formulate the iteration depicted diagram- 
matically in Fig. (5). It yields 



-(a;)<5„,o - a-(a;)<5„,i + a,(a;) f dw' Jz^. 

, LU (jJ 



IT] 

(16) 

Here we encounter the propagator tt{uj) = n(^)o^o(^)- is given by the 
propagator 11, restricted to first order in the tunneling, and we start and 
end in an off-diagonal state. Furthermore, since the parallel tunneling line 
carries a frequency uj the energies of the upper and lower lines are effectively 



Resonant Tunneling and Charging Effects, a Path Integral Approach 13 



shifted relative to one another. Notice that due to the restriction to a two 
state problem there are no diagrams contained in XI^^^q o where a tunneling 
line connects the upper and lower propagator. We can express it by the 
first order self-energy cr{uj) = Yl^-^^ q q{u>) , which is the off-diagonal version 
of the expression known from the first order discussion, with the added 
complication of a parallel tunneling line with frequency uj. 

The irreducible self-energy S is obtained from (p{uj) by connecting the 
tunneling line with the appropriate direction to the upper and the lower 
propagator and adding both contributions (see Fig. (5)). We get for n = 0, 1 

= -2iIm fdwJ2 (f^ni'^) , (17) 
J J, 

while I]„^o follows from the rule S„^o + ^n.i = 0. 

Applying our rules for the calculation of the diagrams we find 

7r(a;) = , a{u;) = - / dcu' , "^"^'^ . . (18) 

^ ' u>-Ao-a{io) ^ ' J io'-LO-iT] ^ ^ 

Here and for the following we introduce the notations a{ij) = a~^{uj) + 
a~{Lo), a'^{Lu) = J2r'^ti'^)j ar{uj) = a+(u;) -|- a~{u:) and = J2r'^o ■ 
Notice that a^{u!) = ar{u})f{±{u: — eVr)) where /(w — eVr) = fr{^) = 
l/[exp(/3(ti; — eVr)) + 1] is the Fermi function of reservoir r. 
We solve the integral equations (16) and obtain 

So,i = -So,o = 27ri^ , Si,o = = 27ri^ (19) 
with \± = j duja'^{uj)\7T{uj)\'^ , \ = j dio\'n{u)\^ . (20) 

Inserting these quantities in the kinetic equation (12) we arrive at the sta- 
tionary probabilities Pq* = A_ and * = A+ . Both probabilities are posi- 
tive and normalized A+ -|- A_ = 1. 



4.3. THE CURRENT 

The expression for the current at time t in the junction r can be written as 

Ir{t) = -2ie{ f dt' J2 «r'^(i - t') sm[Vr,i{t) - VrAt')]) > (21) 

where the expectation value is taken with the reduced density matrix dis- 
cussed above, where the larger of the times is i = t/. If we express the 
expectation value in the charge representation we see that the sin-functions 
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in Eq.(21) describe the transfer of charges at times tf and t'. We, therefore, 
have to evaluate the two real-time correlation functions describing charge 
transfer processes at times t and t' 

C>(t,t') = -^(e-^^(*)e^*'(*')) , C<(t,0 =^(e^*'(*'^e-^^W) . (22) 

In the stationary limit the current can be expressed by 

If = -iejdw {a+{u:)C>{u) + a;{u)C<{uj)} . (23) 

It is further convenient to introduce a spectral density for charge excitations 

A{u) = ^^[C<{u)-C>{u)]. (24) 

The correlation functions and spectral density can be evaluated in the 
approximation which we have used before, 

C<>\u) = +(-)27ri^a,(u;)/(+(-)[u;-eK])|7r(u;)|2 (25) 

r 

and A{uj) = a{Lo) |7r(a;)p . (26) 
The current can then be written as 



/ ^E "^'t?;;^"^ -4(a;)[/.,(a.)-/.(a.)] • (27) 



1, ^^TT 



These results satisfy the conservation laws and sum rules. The current 
is conserved, J2r -^r* — 0' vanishes in equilibrium when = 0. The 
spectral density is normalized / dujA{u!) = 1, and the usual relationships 
between the correlation functions and the spectral density in equilibrium 
are reproduced. The classical result can be recovered if we use the lowest 
order approximation in ao for the spectral density A^^\lo) = 6{uj — Aq). 
We conclude with the observation that quantum fluctuations yield energy 
renormalization and broadening effects, which manifest themselves in the 
spectral density via the real and imaginary part of the self-energy cr{io) 
given in Eq.(18). It will be evaluated in the next section. 

5. Results and applications 

5.1. CHARGE FLUCTUATIONS IN THE SINGLE ELECTRON BOX 

In equilibrium when Vr = Vl, the SET-transistor is equivalent to the single 
electron box. The average excess particle number can be obtained from 
(20) and (26) {n) = J du;f{u;)A{co). In the classical limit A^^\uj) = 5{lo - 
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Aq) one obtains {rf^) = /(Aq) where the energy difference Aq = Ec{^ — 
2CgVg) depends on the gate voltage. Thus, {n'^{QG)) shows a step at 
Qq = CgVg = 1/2, which is smeared by temperature. 

At larger values of or lower temperature we have to include the self- 
energy cr{uj) (18) in the spectral density (26). The two limits, T = and 
|c<;| < T, can be analyzed analytically. In the first case, the spectral density 
has the form 

^(a;)^M , , T = (28) 

^ ^ Ao [uj- A(a;)]2 + [7rA(a;)a(a;)]2 ^ ^ 

where 

The spectral density A(u;) has a maximum at the renormalized energy 
difference A, which is obtained from the self-consistent solution of 

A = A(A) . (30) 

It ftirthcr has a broadening of order nAa where a = a{A). For ira <C 1 the 
broadening can be neglected. In this case our results coincide with the RG- 
analysis of Refs. [14, 11], which shows that the leading logarithmic terms 
are included in our diagram series. 

At finite temperatures |a;| < T, we get 

, A aco coth(T^) , , , , 

^ ^ Ao [lo - A]^ + [Tra Lo coth{^)]^ ' i i- y J 

where 

A= ^^^-^ , a = (32) 

l + 2aoln(|^) l + 2aoln(|^) 

The broadening is of order iraT. It can only be neglected if for vra <C 1. 

As a consequence of quantum fluctuations the step of the average charge 
in the electron box near the degeneracy points is washed out. We neglect 
broadening effects (vra <C 1) and assTimc that the energies of the ground 

state and the first excited state near the degeneracy point depend symmet- 
rically on the distance, ^ =F ^- In this case the partition function is Z = 

2e 4T cosh(^) and the average excess charge, (n) = Qg/& — ^gs^^'^-^i 
becomes 



16 



J. KONIG, H. SCHOELLER, G. SCHON, R. FAZIO 



Depending on the temperature we have to choose the appropriate hmiting 
form for A. At finite temperature the slope at Aq = is 5(n)/5Ao|Ao=o = 
— {4r[l + 2ao In (^^)]2}~^. It is reduced compared to the classical result 
— l/4r. This anomalous temperature dependence can serve as an indication 
for coherent higher order tunneling processes. Below we will demonstrate 
that the linear and nonlinear conductance of a transistor shows more pro- 
nounced deviations from the classical result, making it more suitable for 
experimental verification. 



5.2. CONDUCTANCE OSCILLATIONS IN THE SET-TRANSISTOR 
The linear conductance G of a transistor follows from (27) 

-2 aji{uj)aL{uj) 



G 



du 



-A{u;)f'{u;). 



(34) 



Since the derivative of the Fermi function restricts the integration variable 
to the regime |a;| < T we can use the form (31) and obtain 



G = -27r2 ^ 
h a 



uj/T 



a LO coth (^) 



sinh(w/r) - A]2 + [ttqw coth (^)] 2 



(35) 



where a = + oil and we have to use the finite temperature form for A 
and a. The conductance has a maximum at Aq = 0, given by 



h a 



with a broadening given by 



TT 



— arctan 



1 



2a7r 



(36) 



^ / dAo G(Ao) ^ TT^ T[l + 2ao In (|^)] 

— — • » " — ' 1 ^ ' I 



^max 



TT - 2 arctan [^^!^^ 



(37) 



In the regime 2aoln(^^) <C 1 the peak hight Gmax is constant and 7 
proportional to T, which is the classical result for sequential tunneling. For 
2aoln(^^) ~ 1 and ira <C 1, Gmax and 7 contain logarithmic terms in 
the temperature indicating energy rcnormalization effects due to higher- 
order tunneling processes. For T — > the maximum value as well as the 
broadening vanish, Gmax ~ and 7 ~ TlnT. For 2Q;oln(^^) ~ 1 and 
ttq; ~ 1, we have to retain the finite life-time effects contained in (36) and 
(37). 

The most pronounced signature of quantum fluctuations is contained in 
the differential conductance G{V) = ^^gy-^- In this case the finite voltage 
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Figure 6. The differential conductance in the nonlinear response regime as func- 
tion of the gap energy normalized to the transport voltage V, for T = and 
= — 0.05 {Rl ~ Rr ~ Rt), we consider a symmetric bias and chose 

(1) V/Ec = 0.1 , (2) V/Ec = 0.01 , (3) V/Ec = 0.001 . For comparison, (0) shows the 
result for the classical case obtained from lowest order perturbation theory and the mas- 
ter equation. 

eV provides another energy scale and the renormalization and life-time 
effects are probed over a wider energy range even at zero temperature. The 
T = result obtained from Eq. (27) is plotted in Fig. (6). For comparison we 
also show the result obtained in perturbation theory. (We had quoted the 
expression for the current in this limit in the Introduction.) In the classical 
limit the conductance is nonzero only in the range |Ao| < ^ with vertical 
steps at the edges. The result of Fig. (6) displays clearly the renormalization 
effects and, moreover, the finite life-time broadening. It is clear that at finite 
temperature the conductance is washed out further. These effects should 
be observable in an experiment with realistic parameters. 

6. Interactions in the leads 

We will now address the problem of interactions in the leads by modeling 
the leads as 1-dimensional Luttinger liquids. In contrast to Kane and Fisher 
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[31] we assume that the interactions in the island are sufficiently described 
by the charging energy. We derive a path integral description for this case 
and find modifications of the effective action as compared to the metallic 
case. However, the physics associated with the Coulomb blockade remains 
the same. 

The Hamiltonian of the system is still given by Eq. (1), where the Hamil- 
tonians of the leads describe Luttinger liquids. They are given by [26] 



Hl 



1 

2^ 



(38) 



i.e. they contain a sum of spin {j = cr) and charge (j = p) degrees of 
freedom. The parameters gj describe the interaction strength {gj = 2 in 
the noninteracting limit), and Vj are the velocities of spin and charge ex- 
citations. Wc also introduce the bosonic fields (j)s = <pp + s(f)a- and Og = 
{l/2)[6p + s^ct] for spin up (s = 1) and down (,s = — 1) fermions. These 
fields obey the commutation relation [(j)s{x), 9s'{x')] = (ivr)6(a; — x')6s,s'- 

The fermionic field operators ^ can be expressed in terms of the spin 
and charge degrees of freedom [30] 



X,T 



Jmkpx i^[9p+s0^] i[4>p+scf>a] 



(39) 



where po,s = Ns/2L is the electron density for one spin direction and kp 
the Fermi wave-vector. The bosons fields 6 and (j) can be decomposed in 
terms of scalar bosonic fields bj^q, bj ^ [30] 



ej{x-T) = 

^j{x-T) = 



i 



TT 



q^O 



2qL 



TT 



1/2 



2qL 



1/2 



sign(g)e*«^(el''1^^^6]^^ + g-l^l'^^^^j- _ 



(40) 



The tunneling occurs in the junctions at the points a; = and x = d. 
We obtain an effective action, as in the metallic case, by introducing the 
auxiliary field related to the voltage drop across the junction ^(r) and 
performing a gauge transformation which transfers the phase dependence 
into the tunneling matrix elements. After a cumulant expansion we express 
the partition function as a path integral involving the Euclidean action 



[■P C 



/ dr dr' 2a(r - r') cos[<^(t) - <^(r')]. (41) 
Jo 
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It has the same form as for noninter acting electrons. However, the kernel 
a(r) = ^|r|^iV(0)f , "^ti V (42) 

carries an exponent ry, which is given by the parameters of the LTittinger 
liquid T] = \ + ^[qp + + ^ + Y-]- In Fourier space the kernels arc pro- 
portional to w'^'^^. As a result the I-V characteristics of a normal metal- 
Luttinger liquid junction will also show a power law dependence propor- 
tional to sign(F)y~^. 

We point out that the effective action differs from the one obtained in 
Ref. [31] for the model of a Luttinger liquid with two constriction. The 
physics associated with Coulomb blockade is very similar whether we de- 
scribe the central island as a Luttinger liquid or as a metal island. We finally 
point out that a similar action was discussed in Ref. [29] in connection with 
shake-up effects close to the Fermi energy due to tunneling, and in Ref. [32] 
for junctions coupled to Ohmic baths. 
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